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ABSTRACT 

We present a measurement of the volumetric Type la supernova (SN la) rate out to z ~ 1.6 from the Hubble 
Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera 
for Surveys we discovered 29 SNe, of which approximately 20 are SNe la. Twelve of these SNe la are located 
in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive 
the volumetric SN la rate in four broad redshift bins, finding results consistent with previous measurements at 
z > 1 and strengthening the case for a SN la rate that is >0.6 x 10~ 4 h^ yr _1 Mpc~ 3 at z ~ 1 and flattening 
out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to 
rebin and combine these results with other measurements for increased statistics. Finally, we compare the 
assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that 
different assumptions may induce significant systematic differences between measurements. 
Subject headings: Supernovae: general — white dwarfs — cosmology: observations 
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1. INTRODUCTION 

Type la supernovae (SNe la) are of great importance both 
as astrophysical objects and as cosmological distance indica- 
tors. An accurate knowledge of the rate at which they oc- 
cur (as a function of redshift) is essential for understanding 
both of these roles. Astrophysically, SNe la play an impor- 
tant role in galaxy evolution. They are a major source of 
iron (e.g., Matteucci & Greggio 1986; Tsujimoto et al. 1995; 
Thielemann et al. 1996) and inject energy into the interstellar 
medium (e.g., Dekel & Silk 1986; Scannapieco et al. 2006). 
The SN la rate is necessary to include these effects in galaxy 
evolution models, particularly at high redshifts where much 
of the important galaxy evolution occurs. Cosmologically, 
SNe la are the best-tested method for measuring the scale fac- 
tor of the universe as a function of redshift, with hundreds of 
SNe now employed in the precision measurement of cosmo- 
logical parameters (e.g., Hicken et al. 2009; Amanullah et al. 
2010; Sullivan et al. 2011). Despite their widespread use as 
distance indicators, the process that leads to a SN la is still 
not well understood. SNe la are widely accepted to be the 
end result of a carbon-oxygen (CO) white dwarf (WD) near- 
ing the Chandrasekhar mass limit but how they near that limit 
is not known (see Livio 2001, for a review). This leaves open 
the question of whether high-redshift SNe are different from 
low-redshift SNe in a way that affects the inferred distance. 

Measurements of the change in the SN la rate with redshift 
can be used to distinguish between models of how SNe la 
occur. While there are a variety of SN la progenitor mod- 
els, most fall into two classes: the single degenerate scenario 
(SD; Whelan & Iben 1973) and the double degenerate sce- 
nario (DD; Iben & Tutukov 1984; Webbink 1984). In the sin- 
gle degenerate scenario, the WD accretes mass from a red 
giant or main sequence star that overflows its Roche lobe. In 
the double degenerate scenario, the WD merges with a sec- 
ond white dwarf after orbital decay due to the emission of 
gravitational radiation. Crucially, the delay time between the 
initial formation of the system and the SN explosion is gov- 
erned by a different physical mechanism in the different mod- 
els. This allows us to differentiate between models by mea- 
suring the distribution of the delay times for a population of 
SNe. The shape of this delay time distribution (DTD) de- 
pends on the details of the binary star evolution (particularly 
its evolution through one or more common envelope phases) 
and the specific progenitor model. One method for measur- 
ing the DTD is to correlate the cosmic star formation history 
(SFH) with the the cosmic SN la rate as a function of redshift 
(Yungelson & Livio 2000): the rate as a function of cosmic 
time is simply the cosmic SFH convolved with the DTD. 

The volumetric SN la rate has now been measured in many 
different SN surveys designed to detect and measure SNe at 
z < 1 (e.g. Pain et al. 2002; Neill et al. 2007; Dilday et al. 
2010). With the recently revised rates from the IfA Deep 
survey (Rodney & Tonry 2010), most of these z < 1 mea- 
surements have now come into agreement. In contrast, mea- 
surements at z > 1 have been limited to SN searches in 
the GOODS 28 fields (Dahlenetal. 2004; Kuznetsova et al. 
2008; Dahlen et al. 2008) using the Hubble Space Telescope 
(HST) and ultra-deep single-epoch searches in the Subaru 
Deep Field (SDF) from the ground (Poznanski et al. 2007; 
Grauretal. 2011). These studies have yielded discrepant 
results for the DTD. The first z > 1 measurements by 

28 Great Observatories Origins Deep Survey (Giavalisco et al. 2004) 



Dahlen et al. (2004) (and later Dahlen et al. 2008, with an ex- 
panded dataset) showed a rate that peaked at z ~ 1 and de- 
creased in the highest-redshift bin at z > 1.4. From these 
results the best-fit DTD is one tightly confined to 3-4 Gyr 
with very few SNe la having short delay times (Strolger et al. 
2010). The recent results of Graur et al. (201 1) from the SDF 
show a lower rate at z ~ 1, and a higher rate in the highest- 
redshift bin compared with Dahlen et al. (2008). These results 
are consistent with a flat SN rate at 1 < z < 2. They find 
that the DTD is consistent with a power law with the best-fit 
oc t" 1 ' 1 , implying a significant fraction of short delay time 
(< 1 Gyr) SNe. 

Relative to the HST measurements, the SDF measurements 
cover a much larger volume and therefore have the advantage 
of better statistics in the highest-redshift bin, but HST mea- 
surements hold advantages in systematics. A rolling search 
with HST offers multiple observations of each SN and much 
higher resolution than possible from the ground, useful for 
resolving SNe from the cores of their hosts. These factors 
lead to a more robust identification of SNe la relative to the 
SDF searches where a single observation is used for both de- 
tection and photometric typing. In addition, the Dahlen et al. 
(2008) analysis used spectroscopic typing in addition to pho- 
tometric typing, whereas Grauretal. (2011) uses only pho- 
tometric typing. In general, the very different strategies em- 
ployed make HST measurements a good cross-check for the 
SDF measurements and vice versa. Increasing the statistics in 
HST rate measurements can make this cross-check better and 
improve DTD constraints. At the same time, in comparing the 
measurements it is important to carefully consider possible 
systematic differences, particularly as statistical uncertainty 
decreases and systematics come to dominate. 

In this paper, we address these issues by (1) supplementing 
current determinations of the H5T-based z > 1 SN la rate and 
(2) comparing the effect on results of different dust distribu- 
tions assumed in previous analyses. We use observations from 
the HST Cluster Supernova Survey, a survey to discover and 
follow SNe la in very distant clusters (Dawson et al. 2009, PI: 
Perlmutter, GO-10496). The survey encompassed 189 orbits 
with the Advanced Camera for Surveys (ACS) over a period 
of 18 months. The SN selection and SN typing for the all SNe 
in the survey was presented in Barbary et al. (201 1, hereafter 
Bll), where we calculated the cluster SN la rate from the 
survey. Here, we use a similar methodology to B 1 1 but focus 
on the SNe discovered in the cluster foregrounds and back- 
grounds. The remainder of the paper is organized as follows: 
In §2 we summarize the HST Cluster Supernova Survey and 
the SN discoveries. In §3 we describe the Monte Carlo simu- 
lation used to calculate a rate based on the SN discoveries. In 
§4 we present results and characterize systematic uncertain- 
ties. Finally, in §5 we compare our results to published mea- 
surements. Throughout the paper we use a cosmology with 
Hq = 70 km s" 1 Mpc~\ Q M = 0.3, Ov = 0.7. Magnitudes 
are in the Vega system. 

This paper is one of a series of ten papers (Melbourne et al. 
2007; Barbary et al. 2009; Dawson et al. 2009; 
Morokuma et al. 2010; Suzuki et al. 2011; Ripoche et al. 
2011; Meyers et al. 2011; Hsiao et al. 2011, Bll; This 
work) that report supernova results from the HST Cluster 
Supernova Survey. The survey strategy and SN discoveries 
are described in Dawson et al. (2009), while spectroscopic 
follow-up observations for SN candidates are presented in 
Morokuma et al. (2010). A separate series of papers, ten to 
date, reports on cluster studies from the survey: Hilton et al. 
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Table 1 

Non-Cluster Supernova Discoveries 



ID Nickname z Cluster z Type Confidence 



SNe: Not in Clusters 









1 37 


cc 


plausible 


SN SCP05N10 a 


Tobias 


0.203 


1.03 


CC 


plausible 


SN SCP06C7 




0.61 


0.97 


CC 


probable 


SN SCP06Z5 


Adrian 


0.623 


1.39 


la 


secure b 


SN SCP06B3 C 


Isabella 


0.743 


l.lz 


ll 


probable 


SN SCP06F8 


Ayako 


0.789 


1.11 


CC 


probable 


SN SCP05P9 


Lauren 


0.821 


1.1 


la 


secure^ 


SN SCP06H3 


Elizabeth 


0.85 


1.24 


la 


secure b 


SN SCP06U7 


Ingvar 


0.892 


1.04 


CC 


probable 


SN SCP05P1 


Gabe 


0.926 


1.1 


la 


probable 


SN SCP06G3 


Brian 


0.962 


1.26 


la 


plausible 


SN SCP06C0 


Noa 


1.092 


0.97 


la 


secure 


SN SCP06N33 


Naima 


1.188 


1.03 


la 


probable 


SN SCP06F6 




1.189 


1.11 


non-la 


secure 


SN SCP06A4 


Aki 


1.193 


1.46 


la 


probable 


SN SCP05D6 


Maggie 


1.314 


1.02 


la 


secure 


SN SCP06G4 


Shaya 


1.35 


1.26 


la 


secure b 


SN SCP06X26 


Joe 


1.44? d 


1.10 


la 


plausible 


SNe: Cluster Membership Uncertain 








SN SCP06E12 


Ashley 




1.03 


la 


plausible 


SN SCP06N32 






1.03 


CC 


plausible 



Uncertain to be SN 
SNSCP06M50 C ■•• ••• 0.90 



a Excluded from this analysis due to being inconsistent with an SN la peaking 
< 10 rest-frame days before the first observation. These SNe are excluded 
due to the difficulty of typing SNe found far on the decline. 
b Spectroscopically confirmed. Spectroscopy for SNe SCP05P9, SCP06H3 
and SCP06G4 is reported in Morokuma et al. (2010) 

c Excluded from this analysis due to being within 20" of cluster center. These 

regions are excluded to reduce complications of lensing. 

d Based on a marginal emission line at 9100 A (see Morokuma et al. 2010). 

(2007); Eisenhardt et al. (2008); Jee et al. (2009); Hilton et al. 
(2009); Huang et al. (2009); Rosati et al. (2009); Santos et al. 
(2009); Strazzullo et al. (2010); Brodwinetal. (2011); 
Jee etal. (2011). 

2. SURVEY AND SUPERNOVA DISCOVERIES 

The details of the HST Cluster SN Survey are described in 
Dawson et al. (2009). Briefly, the survey targeted 25 massive 
galaxy clusters in a rolling SN search between July 2005 and 
December 2006. Clusters were selected from X-ray, optical 
and IR surveys and cover the redshift range 0.9 < z < 1.46. 
During the survey, each cluster was observed once every 20 
to 26 days during its HST visibility window (typically four to 
seven months) using ACS. Each visit consisted of four expo- 
sures in the F850LP filter (hereafter zsso)- Most visits also 
included a fifth exposure in the F775W filter (hereafter 1775). 

The process of selecting SN la candidates for the rates anal- 
ysis is described in detail in Bll. We briefly summarize it 
here. Candidates were detected in subtractions of zgso im- 
ages and examined by eye to eliminate obviously false detec- 
tions such as stellar diffraction spikes. A total of 86 candi- 
dates were selected in this phase. Detailed information on all 
86 of these candidates is available from the survey website 29 . 
We generated a full light curve in z§50 an d ^775 for each of 
these candidates and then imposed automated requirements 
on the light curve. These included a requirement on the flux 
in ^775 and the rapidity of the rise and fall of the light curve. 

29 http://supernova.lbl.gov/2009ClusterSurvey/ 



After this step, 60 SN candidates remained. (The selections 
up to this point are accounted for in our calculation of detec- 
tion efficiency.) The remaining candidates were then divided 
into image artifacts (14), AGN (17), and supernovae (29) on 
the basis of the light curve shape and evidence from image 
subtractions. For example, candidates on the cores of bright 
galaxies and showing adjoining positive and negative regions 
in image subtractions are likely to be the result of image mis- 
alignment. With corroborating evidence from the light curves 
of such candidates, they are confidently dismissed as image 
artifacts. Similarly, candidates deemed to be AGN were lo- 
cated on the cores of galaxies and exhibited light curves that 
look nothing like SNe light curves: most rose or fell over peri- 
ods of 100+ days. In general, the continuous light curve mon- 
itoring in the survey made possible to separate these artifacts 
and AGN from SNe with high confidence. For the remaining 
29 candidates deemed to be SNe, we determined a SN type 
and confidence for each. In Table 1 we list the SNe along 
with their host galaxy redshifts, SN types and confidence. We 
omit the eight SNe whose hosts are spectroscopically con- 
firmed cluster members. See Figure 4 of B 1 1 for images, light 
curves and light curve fits of all candidates. A complete de- 
scription of the SN coordinates, typing and confidence level 
(plausible, probable, or secure) is given in Bl 1. Briefly, a se- 
cure SN la has either spectroscopic confirmation, or evidence 
from two sources (early-type host galaxy and light curve) rul- 
ing out other types. A probable SN la is slightly less certain 
than a secure SN la, but still a high-confidence SN la: the 
light curve rules out all core collapse subtypes. A plausible 
SN la has a light curve that is more indicative of a SN la than 
SN CC, but is not sufficient to rule out all core-collapse types. 

In this analysis we use two additional selections not used 
in Bll: (1) First, we eliminate candidates that could only 
be consistent with a SN la if it peaked prior to 10 rest-frame 
days before the first observation. We found that lower-redshift 
(z < 0.9) SNe were detectable even when peaking well before 
the first observation, but that such SNe were extremely diffi- 
cult to type as they were observed only far into the light curve 
decline. We found it most "fair" to eliminate such candidates 
entirely. We include the same selection in our efficiency sim- 
ulations below. Only SCP06L21 and SCP05N10 are rejected 
based on this selection, but both are below the redshift range 
of greatest interest (z > 0.6) and at least SCP05N10 is in- 
compatible with an SN la light curve anyway. This was not 
an issue for Bll because SNe of interest (at z > 0.9) are not 
detectable very far after peak. 

(2) Second, we exclude regions within 20" of cluster cen- 
ters, in order to avoid the most strongly lensed areas in the 
volume behind the clusters. This region is only ^3% of the 
observed field of each cluster. Note that we were careful to 
choose this radius before looking at the radii of any of the 
candidates, in order to avoid biasing ourselves by adjusting 
the radius to conveniently exclude or include candidates. Two 
candidates were excluded as a result: SNe SCP06B3 (16.8" 
from the cluster center) and SCP06M50 (19.4" from the clus- 
ter center). As it happens, these candidates are unlikely to 
be SNe la. SN SCP06B3 is a "probable" SN CC, while SN 
SCP06M50 is possibly not a SN at all and may be hosted by 
a cluster member galaxy, making its position near the cluster 
center unsurprising. The exclusion of this region is taken into 
account in our simulations (§3). The effect of lensing on the 
remaining portions of the fields are discussed in §4.2. 

The systematic uncertainties associated with the determina- 
tion of SN type and redshift for the remaining candidates are 
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K09 Dust 
Dust Model A 
Dust Model B 
Dust Model C 
Minimal Dust 




Figure 1. Left Panel: Stretch distribution used for simulated SNe (solid black line) and the stretch distribution of first-year SNLS z < 0.6 SNe (grey histogram) 
from Astier et al. (2006). Right Panel: Color distribution used for simulated SNe (solid black line), based on the K09 distribution of host-galaxy extinction. The 
grey histogram shows the color distribution of the first-year SNLS z < 0.6 SNe. The other four lines show alternative color distributions used to assess the 
possible systematic uncertainty due to different distributions of host galaxy dust extinction (see §4.3). 



addressed in 84. 1 . 



3. RATE CALCULATION 

We calculate the SN la rate in redshift bins using what has 
become a standard method in rate calculations: The number 
of SNe la per unit time per comoving volume is estimated in 
the redshift bin z\ < z < z 2 by 



TZ(zi < z < z 2 ) = 



NsN IaCgl < Z < Zi) 



(i) 



where A*sn \ & {z\ < z < z 2 ) is the number of SNe la dis- 
covered between redshifts z\ and z%, and the denominator is 
the total effective time-volume for which the survey is sensi- 
tive to SNe la in the redshift range z\ < z < z<i. T(z) is 
the effective visibility time (also known as the "control time") 
and is calculated by integrating the probability of detecting a 
SN la as a function of time over the active time of the survey. 
T{z) depends on the dates and depths of observations, as well 
as the specific requirements for selecting SNe. The factor of 
1/(1 + z) converts from observer-frame time to rest-frame 
time at redshift z. The last two terms in the denominator rep- 
resent the volume comoving element between z and z + dz 
observed in the survey. * s me comoving volume of a 

spherical shell of width dz. is the solid angle observed in 
the survey, in units of steradians. (Q/Att is the fraction of the 
spherical shell we have observed.) Finally, the average red- 
shift of the bin, weighted by the volume effectively observed, 
is given by 



■C*r(*)ik%(*)<fe 



(2) 



As in Bl 1, we use an effective visibility time that depends 
on position, as observation dates and depths vary within each 
observed field. That is, in Equation (1) we make the substitu- 
tion 



T(z)e 



T(x,y,z)dxdy. 



(3) 



■'■■;/ 



T(x, y, z) is calculated by simulating SN la light curves at 
different positions, redshifts and times during the survey, and 
determining the probability that each simulated SN would be 



detected and counted in our SN sample. We pass each simu- 
lated SN through the same automated selections used to select 
the 60 candidates in our initial sample. Additionally, we dis- 
card simulated SNe peaking prior to 10 rest-frame days before 
the first observation, as discussed in the previous section. 

We characterize the diversity of SN la light curves as a two- 
parameter family (stretch s and color c) with an additional 
intrinsic dispersion in luminosity. The absolute magnitude of 
each simulated SN is set to 



M B = -19.31 - a(s -l)+f3c + I 



(4) 



where —19.31 is the magnitude of an s = 1, c = SN la 
in our assumed cosmology (Astier et al. 2006), a = 1.24, 
(3 = 2.28 (Kowalski et al. 2008), and / is an added "intrin- 
sic dispersion", randomly drawn from a Gaussian distribution 
centered at zero with a = 0.15 mag, as seen in Kowalski et al. 
(2008). This produces a set of simulated SNe that closely 
matches the distribution seen in Kowalski et al. (2008). To 
calculate the flux of each simulated SN in the observed z S50 
and 1775 filters, we use the Hsiao et al. (2007) spectral time 
series template. 

The main difference from Bll is that we use distributions 
for stretch and color that are representative of SNe in the field 
rather than in clusters. For stretch, we look to the observed 
stretch distribution of the first-year sample from the Super- 
nova Legacy Survey (SNLS; Astier et al. 2006), cut at z < 0.6 
to limit Malmquist bias (Figure 1, grey histogram). We as- 
sume this sample is complete, fit a smooth curve to the dis- 
tribution (same panel, solid line), and use this for simulated 
SNe. 

For color, we cannot assume the SNLS sample (Figure 1, 
histogram in right panel) is complete even at z < 0.6, as 
highly reddened SNe will have been missed. The standard 
picture today is that the observed distribution of SN colors is 
due to a combination of both intrinsic SN color variation and 
host galaxy extinction (Guy et al. 2010; Chotard et al. 201 1). 
Both of these are expected to introduce a color that corre- 
lates with SN luminosity, possibly with different strengths (f3 
is typically found to be smaller than the canonical value of 
Rb = 4.1 for Milky Way dust). In order to capture both ef- 
fects with a single color distribution and a single (3, we work 
backwards from the desired host galaxy extinction distribu- 
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tion. We wish to achieve a host galaxy extinction distribu- 
tion of P(A V ) oc cxp(-^4y/0.33) with A v > 0, the best-fit 
value for host-galaxy SN extinction in the SDSS-II SN Sur- 
vey (Kessler et al. 2009, hereafter K09). To do this we use a 
color distribution of P(c) oc exp(— (/3 — l)c/0.33), because 
A v is related to c via A v = Rv x E(B — V) w (J3 — 1) x c. 
This ensures the desired Ay distribution is obtained for any 
given value of (3. The full P(c) distribution is then a convo- 
lution of this P(c) distribution from host galaxy dust and the 
intrinsic distribution of SN color (assumed to be Gaussian). 
The Gaussian parameters of the intrinsic distribution are cho- 
sen so that the full convolved P(c) distribution matches the 
observed SNLS c distribution at c < 0.3. The result is a 
P{c) distribution (black line in right panel of Figure 1) that 
matches the SNLS sample where we expect it to be complete 
(c < 0.3) and also has the desired behavior at large extinction 
based on our assumed P(Ay) distribution. We use this distri- 
bution in our simulation. In §4.3 we assess the systematic un- 
certainty associated with host galaxy dust by using alternate 
color distributions obtained using the same method, but dif- 
ferent P{Ay) distributions (other curves in the same panel). 

T(x, y, z) is calculated in bins of 100 x 100 pixels (5" x 
5") in position. In each 5" x 5" region, we simulate 50 
SN light curves with random parameters and random position 
(within the bin) and take the average effective visibility time 
of the 50 SNe (-80,000 SNe per field). Summing over all 
areas observed in all 25 fields yields T(z). In doing so we 
exclude regions within 20" of cluster centers, as discussed in 
the previous section. We calculate T(z) at intervals of Az = 
0.05 inredshift. 

4. RESULTS AND SYSTEMATIC UNCERTAINTY 

Figure 2 (top panel, black line) shows the product of the 
observer-frame effective visibility time and the area (T(z)0 
from Eq. 1) as a function of SN redshift. For reference, 
the horizontal dotted line shows an approximate calcula- 
tion of this value, multiplying the area of the ACS field 
(11.65 arcmin 2 ) by the time difference between 10 days be- 
fore the first observation and 10 days after the last observa- 
tion. In reality the area actually observed is slightly more 
complicated and SNe are detected over a slightly larger time 
range. From z = 0, the effective visibility time actually in- 
creases slightly out to z — 0.5 as SN light curves are time- 
dilated and are thus visible for longer. Afterwards, we begin 
to miss SNe that peak during the observations. In the lower 
panel of Figure 2, we convert to the rest-frame time-volume 
observed in each redshift bin of Az = 0.05 using the assumed 
cosmology. 

Table 2 shows the results, in bins of Az = 0.4 (comparable 
to Dahlen et al. 2008) and also in bins of Az = 0.5 (compara- 
ble to Graur et al. 201 1). The numerator of Equation (1) (third 
column in Table 2) is the number of SNe la described in §2. 
The denominator of Equation ( 1 ) (fourth column in Table 2) is 
obtained by summing Figure 2 (lower panel) over the redshift 
bin of interest. We now discuss the systematic uncertainties 
associated with lensing, SN type determination, host-galaxy 
dust, SN properties, and galaxy number density variations. 

4. 1 . Type determination 

The uncertainty in SN type in the survey is quite small, 
thanks to the cadenced nature of the survey and excellent 
spectroscopic follow-up. Consider the candidates designated 
as SN la: all three SNe la at z < 0.9 are spectroscopically 
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Figure 2. Top Panel: The observer-frame effective visibility time multiplied 
by observed area, as a function of supernova redshift. The horizontal dot- 
ted line shows the area of the ACS field multiplied by the time spanned by 
the observations in each cluster. Bottom Panel: The rest-frame volume-time 
searched in each redshift bin of Az = 0.05. In each panel, the black line 
shows our main result for the effective visibility time, based on simulations 
using the K09 dust distribution. The green, red, blue and cyan lines show the 
results for alternative dust distributions. [Data in this figure is available as a 
machine-readable table.] 



confirmed. At z > 0.9, any SN bright enough to be detected 
is overwhelmingly likely to be Type la due to the faintness of 
core-collapse SNe relative to SNe la (e.g., Dahlen et al. 2004; 
Li et al. 201 1; Meyers et al. 201 1). Furthermore, while "prob- 
able" candidates are not as certain as "secure" candidates, this 
is still a fairly high-confidence type determination: A "prob- 
able" SN la means that a SN la light curve template has a 
X 2 P-value that is 10 3 times larger than any SN CC value. 
A Bayesian analysis would therefore yield a type uncertainty 
close to zero for such candidates, regardless of the prior used. 

The "plausible" candidates are perhaps the only candi- 
dates with significant type uncertainty. It is difficult to pre- 
cisely quantify this uncertainty. Instead, we provide conser- 
vative bounds in a manner similar to Dahlen et al. (2008) and 
Sharon et al. (2010): we first assign a lower limit to the num- 
ber of SNe la discoveries by assuming that all "plausible" 
SNe la are in fact SNe CC. We then assign an upper limit 
by assuming that all "plausible" SNe CC are in fact SNe la. 
These limits are shown in Table 2 as the second confidence in- 
terval for Nsn ia- The corresponding systematic uncertainty 
in the SN rate is shown in Table 3. 

For the two candidates without spectroscopic host redshifts, 
we assign a redshift range consistent with the SN light curve 
and/or host galaxy photometry, as follows: For SCP06E12, 
we use the range 0.8 < z < 1.2. As there is uncertainty about 
both the type and cluster membership, we count SCP06E12 as 
0.5 ± 0.5 field SNe Ia. The situation is similar for SCP06N32: 
the light curve is consistent with an SN Ibc at z — 0.9, but also 
with an SN Ia at z — 1.3. We therefore assign a redshift range 
of 1.1 < z < 1.5 and count it as 0.5 ±0.5 field SNe Ia. These 
two SNe are assigned to redshift bins in Table 2 according to 
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Table 2 

Results: SN la Rate 



Redshift bin 


2 




Denom b 


Rate c 








bin width A 2 = 


0.4 




0.2 < z < 0.6 
0.6 < z < 1.0 
1.0 < z < 1.4 
1.4 < z < 1.8 


0.442 
0.807 
1.187 
1.535 


n m +1.16+0.00 

u - uu -o.oo-o.oo 

c 9 c+2.69+0.25 
°' ZO -2.00-1.25 
c K o+2.77+0.63 
O.Oo_ 2 .08-0 .63 
i i q+1.56+0.12 
± - lz -0.79-1.12 


2.332 
4.464 
4.243 
1.453 


n nn + 0.50 + 0.00 
UiUU -0.00-0.00 
i ls +0.60+0.44 

^"-o. 45-0. 28 

-, OO+0.65+0.69 
l-oo_ .4g_o .26 
n 77 +1.07+0.44 
u -' '-0.54-0.77 






bin width A 2 = 


0.5 





0.5 < 2 < 1.0 
1.0 < 2 < 1.5 
1.5 < 2 < 2.0 



0.766 
1.222 
1.639 



U,UU -0.00-0.00 
c ocr+2.69+0.25 
d - zd -2.00-1.25 
f, 7 c-+2.99+0.75 
O- '0_2.3i_i.75 
n m +1.16+0.00 
u,uu -0.00-0.00 



1.624 
5.321 
4.906 
0.890 



n nn+ u - ' 1 + u ' uu 

u ' uu -0.00-0.00 
n QQ+0.51+0.33 
u -aw_0.38-0.24 
i OO+0.61+0.71 
1-O0_0.47-0.43 
fl. 30+0. 00 



0.00" 



o.oo-o.oo 



a Number of SNe la in bin. The first and second confidence intervals repre- 
sent the Poisson uncertainty and the uncertainty in type determination, respec- 
tively. The non-integer number of SNe in each bin is attributable to the two 
candidates without spectroscopic redshifts. These candidates are assigned 
redshift ranges that are spread over multiple bins. 

b Denominator of Equation (1): the total rest frame time-volume searched in 
n4 h~ 3 Mnr-3 



The rate in units of 10 4 h% n yr 1 Mpc 3 . The first and second confidence 



this bin, having units 10 h 7Q yr Mpc 

^o y r 

intervals represent the statistical and systematic uncertainty, respectively. The 
statistical uncertainty is entirely due to Poisson uncertainty in Ns~n la- The 
systematic uncertainty is the typing uncertainty in A^n la and systematic 
uncertainties in "Denom" (described in text) added in quadrature. 



Table 3 

SN Rate Uncertainties in Percentage 





statistical 




systematic 




Redshift bin 


Poisson 


Typing 


Dust 


Luminosity 


0.6 < 2 < 1.0 
1.0 < 2 < 1.4 
1.4 < 2 < 1.8 


+51/-38 
+49/- 36 
+138/-69 


+4/-23 
+11/-11 
+11/- 100 


+37/- 3 
+49/- 14 
+39/- 19 


+2/-2 
+10/-7 
+38/-23 



Note. — Percentages are not reported for the 0.2 < z < 0.6 bin because 
there were no SNe detected in this bin. 

the degree of overlap between the redshift range of the SN and 
the range of the redshift bin. 

Finally, note that we have classified the highest-redshift SN 
SCP06X26 as a lower-confidence "plausible" SN la, despite 
the fact that any SN detected at z = 1.44 is overwhelmingly 
likely to be Type la. This conservative approach was taken 
because the spectroscopic host galaxy redshift of SCP06X26 
is based on a single low signal-to-noise emission line; as a 
result, the redshift and type are less certain than they would 
otherwise be. The "plausible" designation therefore includes 
the possibility that the SN is at a lower redshift. As a result, 
the confidence interval on the number of detected SNe la in 
the 1.4 < z < 1.8 bin is [0, 1]. Still, the light curve of 
SCP06X26 is consistent with a typical SN la at z = 1.44, so 
this remains the best estimate. 

4.2. Lensing due to Clusters 

The presence of a massive galaxy cluster in each of the 
25 observed fields presents a complication for measuring 
the volumetric field rate. A cluster will preferentially mag- 
nify sources behind it (increasing the discovery efficiency of 
SNe), and will also shrink the source plane area behind 
the cluster (decreasing the number of SNe discovered) (e.g., 



Sullivan et al. 2000; Goobar et al. 2009). Fortunately the ef- 
fect of lensing on the calculated rates in this survey is small for 
three reasons. First, the high redshifts of the clusters means 
that the volume of interest in the cluster backgrounds is close 
to the clusters and therefore not lensed very efficiently. Sec- 
ond, we have already excluded from the analysis the central 
20" of each field, where lensing effects are the largest. Third, 
at any given radius the two effects (magnification and source 
plane area shrinkage) are opposing in terms of number of SNe 
discovered. 

We have calculated the magnitude of each lensing effect on 
the remaining outer regions using a simple lensing model: We 
assume each cluster has a mass of M200 = 4 x 10 14 M 
(the approximate average mass in our sample, as reported 
by Jee et al. (2011) and an NFW (Navarro et al. 1997) mass 
profile. We distribute clusters according to their redshifts 
and calculate the lensing effect on the 25 annular regions 
20" < r < 100" around the clusters. The distribution of 
magnification in these regions as a function of source redshift 
is shown in Figure 3. The magnification is quite small: even 
at a source redshift of z = 1.8, most of the area is magni- 
fied by less than 10%. As a rough estimate of the effect on 
the derived rates, we show the average magnification m avg 
for each source redshift, and the effect such a magnification 
would have on the detectability of SNe at this redshift. To 
calculate the effect on the detectability, we increase the lu- 
minosity of all SNe in our Monte Carlo efficiency simulation 
by a factor m avg and recalculate the denominator of Equa- 
tion (1). More luminous SNe causes "Denom" to increase, 
corresponding to a decrease in the inferred SN rate. The ef- 
fect is only a few percent at 2 < 1.4, where the survey is most 
sensitive, but starts to increase steeply towards z = 1.8. In 
Figure 4 we show the decrease in the source area as a func- 
tion of redshift, which translates directly into a decrease in the 
effective visibility time x area. The decrease is nearly linear 
with redshift past z = 1.2, reaching ~13% at z = 1.8. 

We conclude from these simulations that the two effects 
cancel to within a few percent of the total rate, over the red- 
shift range of interest: At z = 1.4, magnification increases SN 
detectability by ^3% and source-area reduction decreases de- 
tectability by ^6%. At z = 1.6, the increase is ^12%, and 
the decrease is ~40%. At z = 1.8 the increase overwhelms 
the decrease (~27% versus ~43%), but there will be very few 
SNe detected beyond z ~ 1.6 (see Figure 2). Therefore, we 
have not made an adjustment for these effects. Furthermore, 
the size of each effect is much smaller than other sources of 
systematic uncertainty considered below. For example, the 
average magnification at z = 1.8 is only ~4.08 (—0.08 mag), 
whereas below we consider the effect of changing the lumi- 
nosity of all SNe in our simulation by ±0.2 magnitudes. So, 
we do not assign a specific systematic uncertainty to lensing 
effects. 

4.3. Dust Extinction 

The degree to which SNe are affected by host galaxy dust 
extinction is perhaps the largest systematic uncertainty in 
SN la rate studies. Here, we consider alternatives to the ex- 
tinction distribution used in §3 and evaluate the effect on the 
results. Specifically, we reproduce three P(Ay) distributions 
considered in Dahlen et al. (2008) and Neill et al. (2006). In 
each of these models, as for our main result, we constrain 
P(A V ) = Ofor A v < 0. 

The first, "Model A," is used for the main result in 
Dahlen et al. (2008). It is based on the model of Hatano et al. 
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Figure 3. Results of our lensing simulation: the magnification distribution in 
regions at radius 20" < r < 100" in all 25 cluster fields (the approximate 
extent of the regions used in the rate analysis). The four panels correspond 
to source redshifts of z s = 1.2, 1.4, 1.6, 1.8. For each source redshift, the 
average magnification m avg is given. Under the approximation that all SNe 
are magnified by m aV g, A(Dcnom) shows the corresponding increase in 
the denominator of Equation (1). Note that this change in Denom considers 
only magnification; the effect of source-plane area shrinkage is considered 
separately in Figure 4. At z < 1.4, where the survey's sensitivity is greatest, 
magnification from lensing has only a < 3% effect on the detectability of 
SNe. 
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Figure 4. True source-plane area relative to the observed area in our lensing 
simulation, as a function of source redshift. The relative area is for regions 
at radius 20" < r < 100" in each of the 25 cluster fields (the approximate 
extent of the regions used in the analysis). The relative source area corre- 
sponds directly to the change in the denominator of Equation (1). The effect 
is opposite in sign to the effect in Figure 3. 



(1998), constructed to estimate extinction in local disk galax- 
ies. The distribution in Figure 3 of Dahlen et al. (2008) is well 
approximated by 



P(A V ) = 0.61- 



-A v /2 



0.39- 



-A v /0.07 
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Figure 5. Host galaxy dust extinction distributions, illustrating the behavior 
of the distributions at large Ay. The K09 distribution is used for our main 
result. Models A, B, and C are similar to the models of the same name ex- 
amined in Dahlen et al. (2008) and are based on results from Hatano et al. 
(1998), Riello & Patat (2005) and Neill et al. (2006), respectively. These al- 
ternative distributions are used here to investigate possible systematic uncer- 
tainty due to host galaxy dust. This figure can be compared to Figure 3 of 
Dahlen et al. (2008). 

This distribution is shown in Figure 5. The second distribu- 
tion we consider, "Model B," is used by Dahlen et al. (2008) 
as an alternative distribution. It is based on the models in 
Riello & Patat (2005), which are aimed at generalizing the 
Hatano et al. (1998) models to a variety of dust properties and 
a variety of spiral Hubble types. Here we approximate the dis- 
tribution by 



P{A V ) = 0.35S(A V ) +0.40 



0.25e~ 



-A^/(2x0.6 2 ) 



0.6V27T 



(6) 



where 5 is the Dirac delta function, used here to assign 35% 
of SNe to the lowest host galaxy extinction bin. The third 
distribution we consider, labeled "Model C," is used in the 
rate analysis of Neill et al. (2006). It is given by 



-A^/(2x0.62 2 ) 



0.62a 



(7) 



(5) 



All three of these distributions are reproduced in Figure 5, 
and the corresponding distributions of SN color are shown in 
Figure 1 (right panel). In addition to these three distributions, 
we also consider a distribution with minimal dust, where we 
assume the SNLS z < 0.6 sample is complete and fit it with a 
skewed Gaussian distribution. The effective visibility time for 
each dust model is shown in Figure 2 and the corresponding 
SN rate results are shown in Figure 6 (left panel). 

Of all the models, Model A produces the most strikingly 
different results for the effective visibility time. Even in the 
lowest redshift bin (0.2 < z < 0.6) it implies that 10% of SNe 
are missed due to dust, relative to the K09 model. In the 0.6 < 
z < 1.0 redshift bin, it yields effective visibility times lower 
by 27%, while Models B, C and the minimal dust model result 
in changes of only —7%, —3%, and +3% respectively. This is 
unsurprising: In Model A, 26% of SNe have host galaxy ex- 
tinctions Ay > 2 while this fraction is < 4% in Model B and 
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Figure 6. The volumetric SN la rate in four redshift bins (points with error bars) of width Az = 0.4. The error bars represent the statistical-only uncertainty. 
The black line shows the rate calculated in a moving bin of width Az = 0.4 (shaded grey regions represent uncertainty). Note that the points with error bars are 
uncorrelated errors (using non-overlapping bins), while the uncertainty in the moving bin is correlated from point to point. Left Panel: The green, red, blue and 
cyan lines show the rate (with no uncertainty) assuming alternative SN color distributions. Right Panel: The red and blue lines show the rate assuming that all 
SNe are brighter or dimmer by ±0.2 mag. 



< 1% in models K09 and C. In the 1.0 < z < 1.4 bin Model 
A has the largest effect: —33% compared to the K09 model. 
As the work of Riello & Patat (2005, Model B) is aimed at 
generalizing Hatano et al. (1998, Model A) for use at higher 
redshifts, Model B may be viewed as more applicable. Along 
these lines, Cappellaro et al. (1999) have noted that using the 
Hatano et al. (1998) model appears to over-correct SN rates. 

Still, for the systematic uncertainty associated with the 
choice of dust model we take a conservative approach, us- 
ing the full range of these models. We use the minimal dust 
model to obtain the lower limit on the rate and Model A for 
the upper limit. This confidence range is shown in Table 3. 

4.4. More Dust at High Redshift? 

Several studies have pointed out that extinction is likely 
to increase with redshift, due to increasing star formation in 
dusty environments (Mannucci et al. 2007; Holwerda 2008). 
The potential effect on the SN la extinction distribution is 
difficult to quantify as it depends not only on the relation 
between star formation and dust ejection but also the SN la 
delay time distribution. However, we estimate that any pos- 
sible redshift-dependence is well-encompassed by the con- 
servatively large range of extinction models we use above. 
Rowan-Robinson (2003) estimated that the average extinction 
(Ay)(z) peaks at z ~ 1, at a value ^0.15 mag higher than lo- 
cally. The effect of such an additional dimming on our rates at 
z ~ 1 is approximately 10%, whereas the extinction distribu- 
tion uncertainty already included above is ^50%. Similarly, 
Grauret al. (2011) estimated the fraction of missing SNe to 
be 5-13% in the redshift range 1 < z < 2 based on the work 
of Mannucci et al. (2007). As the difference induced by our 
use of the K09 extinction model rather than model C used by 
Graur et al. is already much greater than this, we do not make 
an explicit correction for increasing dust at high redshift. 

4.5. Other SN properties 

In addition to the choice of the host galaxy dust extinction 
distribution, other assumptions about SN properties can af- 
fect the derived rates. These include the absolute magnitude, 



the stretch distribution, and the spectral time series template. 
Fortunately these properties are well constrained. For exam- 
ple, shifting the stretch distribution by As = ±0.05 would 
be inconsistent with the observed distribution, and similarly 
for the color distribution. The average absolute magnitude 
of SNe la (in our cosmology) is constrained to much better 
than 0.1 mag (Amanullahet al. 2010). To first order, chang- 
ing any of these assumptions affects the derived rate by in- 
creasing or decreasing the luminosity of the simulated SN and 
thereby increasing or decreasing the effective visibility time. 
Therefore, to estimate the effect of such changes on the re- 
sult, we simply shift the absolute magnitude of the simulated 
SNe la. A shift in stretch of As = ±0.05 is equivalent to a 
magnitude shift of AMg = qAs = 0.08. Similarly, a shift 
in color of Ac = ±0.05 is equivalent to a magnitude shift 
of AMb = /3Ac = 0.11. Conservatively then, we use a 
±0.2 mag shift to jointly capture these uncertainties and un- 
certainty in the absolute magnitude. The effect on the results 
is shown in Figure 6 (right panel) and is generally comparable 
to or smaller than the effect from the extinction distribution 
uncertainty (see Table 3). 

4.6. Cosmic Variance and galaxy-cluster correlation 

There are various effects that could increase or decrease the 
density of galaxies in the observed fields relative to the cosmic 
average, potentially biasing the volumetric rates. We have 
estimated all such effects to be small (< 5%). We briefly 
discuss three of these effects. 

(1) Masking volume surrounding clusters. Any SN occur- 
ring in the volume around a cluster (within |Az| < 0.015) 
but not associated with the cluster would be mistakenly as- 
signed to the cluster. This volume is therefore effectively not 
counted for this analysis. We estimate that this decreases the 
total volume in the 0.6 < z < 1. redshift bin by ss 1.5% (due 
to the presence of 5 clusters) and in the 1.0 < z < 1.4 bin by 
«6% (due to the presence of 20 clusters). Because the next 
effect compensates somewhat for this effect, we do not make 
an explicit correction in the result. 

(2) Cluster-galaxy correlation. There may be an increase 
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Figure 7. Volumetric SN la rates from the HST Cluster Supernova Survey (red points) compared to key rates from the literature. For measurements with 
two error bars, including ours, the inner and outer error bars represent the statistical (Poisson) and total (statistical + systematic) uncertainties, respectively. 
Measurements with a single error bar (Ku08 and Gil) are Bayesian-based analyses where the error bar encompasses both statistical and typing uncertainties. 
Some measurements have been offset slightly in redshift for visibility. Note that we have conservatively reported our measurement in the highest-redshift bin as 
having a lower limit of zero due to redshift uncertainty in the single SN in this bin. 



in the density of field galaxies (relative to the cosmic average) 
in the volume surrounding the cluster due to the presence of 
the cluster itself. However, we estimate that this effect should 
be small outside of the volume discussed in the previous ef- 
fect. At z = 1, a redshift difference of Az = 0.015 corre- 
sponds to a comoving distance of ^36 Mpc. Density fluctu- 
ations on these scales are in the linear regime. At 36 Mpc, 
the correlation function of massive galaxies is £ < 0.2 (e.g., 
Eisenstein et al. 2005; White et al. 201 1). Even if the average 
over-density at scales 36 Mpc < r < 100 Mpc were 0.2 (an 
overestimate), this would only be a 2-3% effect on the total 
galaxy density in a bin of width Az = 0.4. 

(3) Cosmic variance. Cosmic variance could affect the re- 
sults if there is an under-density or over-density of field galax- 
ies in the observed regions (relative to the cosmic average). In 
Bl 1, we estimated the size of this effect in the GOODS fields, 
finding it to contribute a systematic uncertainty of 4-6% based 
on the cosmic variance calculator of Trenti & Stiavelli (2008) 
and the difference between the GOODS-N and GOODS-S 
fields (see also discussion in Dahlen et al. 2008). The effect is 
even smaller for the 25 fields in this study: here we have more 
than two widely separated fields, providing a better sampling 
of cosmic variance. 

5. DISCUSSION 

The results from this study are available as a machine- 
readable table from the Journal or the survey website 30 . We 
include the product of the effective visibility time and ob- 
served area, as a function of redshift (Figure 2), for a variety 

30 http://supernova.lbl.gov/2009ClusterSurvey/ 



of assumptions about SN properties and host galaxy dust dis- 
tributions. With these data and the SN candidate list, the rates 
from this survey can be recomputed in any arbitrary redshift 
bin and for any of these assumptions. This will make it easy to 
combine these results with other measurements for increased 
statistical power. 

5.1. Comparison to Other High-Redshift Measurements 

In Figure 7 we compare our results to an assortment of 
other volumetric SN la rate measurements. (All measure- 
ments have been corrected to our assumed cosmology.) Our 
results are generally consistent with the three published mea- 
surements at z > 1: Kuznetsova et al. (2008, hereafter Ku08), 
Dahlen et al. (2008, hereafter D08), and Grauretal. (2011, 
hereafter Gil). D08 and Gil supplant earlier results from 
Dahlen et al. (2004) and Poznanski et al. (2007), respectively. 
The Ku08 and D08 measurements are based on SN searches 
in the HST GOODS fields, with Ku08 being an indepen- 
dent analysis of a subset of the data used in D08. These SN 
searches used ACS to cover the GOODS fields with a 45 day 
cadence and triggered followup (imaging and spectroscopy) 
of SN candidates. The D08 analysis uses a SN typing method 
based on both spectroscopy and photometry (similar to the ap- 
proach used here) while Ku08 use a photometric-only pseudo- 
Bayesian approach to typing. The Gl 1 measurement is based 
on "single-detection" searches in the Subaru Deep Field. Gil 
also use a pseudo-Bayesian typing approach, but use a single 
detection with observations in three filters, rather than multi- 
ple detections with observation in (typically) two filters as in 
Ku08. 
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Table 4 

Rate Comparison Using Consistent Extinction Distributions 





0.6 < 2 < 1.0 


1.0 <z < 1.4 


1.4 < z < 1.8 


This work using extinction model A 
D08 


^ — 0.62 
J " ou — 0.27 


^ — 0.74 
i o 9 +0.36 
J -° -0.29 


1.08 + ^° 

42+ - 39 
-0.23 


This work using minimal extinction 
model 


1 14 +0.58 
J - L *— 0.43 


1 14 +0.56 
J-- 1 -0.42 


62+ - 86 
u - DZ -0.44 


Ku08 


qq+0-25 


7 t;+0.35 
u "'°-0.30 


12+ - 58 




0.5 < z < 1.0 


1.0 <z < 1.5 


1.5 < z < 2.0 


This work using extinction model C 
and +5 - 13% high- z correction 


1 oi +0 - 52 

1,u± -0.39 


1 69+ - 75 


00+ 1 92 
u - uu -o.oo 


Gil 


7O+0-33 


fl R4+ - 25 
U - 84 -0.28 


1 02+ 054 



Note. — Rate in units of 10 4 h^ yr 1 Mpc 3 . Confidence intervals are statistical (Poisson). 



5.2. r/je Effect of Different Extinction Distributions 

Although our results generally agree with all three previ- 
ous studies (Ku08, D08, Gil), it is interesting to compare in 
more detail the assumptions used in each study. Here we focus 
specifically on differences in the assumed extinction distribu- 
tion, which we find to be the dominant systematic. Unlike 
systematic uncertainties arising from SN typing, systematic 
uncertainty due to assumptions about SN properties will be 
common between rate studies, potentially leading to a system- 
atic offset between results if different assumptions are used in 
different studies. 

For example, our results appear very similar to D08 in the 
two mid-redshift bins (within 10% in both bins). However, 
D08 assume a different extinction distribution in their simula- 
tions than we do. If we assume the same distribution (extinc- 
tion model A), our derived rates in these two bins are 1.61 and 
1.99 xlO" 4 h% , compared to 1.30 and 1.32 xl0~ 4 h% in 
D08. (Conversely, if D08 had used the K09 extinction model, 
their rates would have been lower than ours.) In other words, 
we actually found more SNe la than one would have predicted 
based on D08 (but not inconsistently so given the large Pois- 
son uncertainties). 

Note that D08 also compare Models A, B and C in assess- 
ing their systematic uncertainty but do not find the large dif- 
ferences that we find here. They find that Model B produces 
rates that are < 10% lower than Model A (only ~4% in the 
highest redshift bin), and that Model C has even less of an 
effect, whereas we find that the difference beyond Models A 
and B is approximately 26% in the central two bins. These 
different findings regarding the systematic uncertainty due to 
extinction are surprising given that both studies are based on 
HST searches the same band to similar depths. One possi- 
ble explanation is the difference in cadence (~23 days here 
versus ~45 days in GOODS). We checked the impact of a 
longer cadence by rerunning our simulations ignoring every 
other epoch in our search. We find that the effect of different 
extinction distributions is unchanged to within a few percent. 
That is, the shorter ^23 day cadence is not significantly more 
sensitive to highly reddened SNe than the ~45 day cadence. 
The slight enhancement in reach that a short cadence affords 
is insignificant compared to the long Ay tail of, for exam- 
ple, Model A. It is possible that the difference in estimated 
systematic uncertainty could arise from some other detail of 
the efficiency simulations, but a full resimulation of the D08 
efficiency calculations is beyond the scope of this paper. 

In Table 4 we have recomputed our rates using extinction 



assumptions similar to D08, Ku08 and Gil. Ku08 use only 
two discrete values for Ay. The values, 0.0 and 0.4, are rela- 
tively small, so we compare using our "minimal dust" model. 
Gil use the distribution of N06 (Model C) with a +5-13% 
correction in the redshift range 1.0 < z < 2.0 as discussed 
earlier (§4.4). Table 4 serves two purposes: (1) it aids direct 
comparison between each of these studies and our result, and 
(2) it illustrates how much the assumed extinction distribution 
affects our result. 

In light of the large systematic differences due to dust, it is 
vitally important to use caution when comparing rate results 
from studies that use different dust assumptions. In particu- 
lar, systematic offsets from dust assumptions could affect the 
shape of the derived SN la delay time distribution. The DTD 
shape is obtained not from the SN rate itself, but from the 
change in the rate with redshift (see, for example, Figure 2 of 
Horiuchi & Beacom 2010). To measure the change, one typi- 
cally must compare surveys covering different redshift ranges. 
Comparing low- and high-redshift measurements that correct 
their rates using different extinction distributions will induce 
a systematic error on the slope of the SN rate and thereby the 
DTD shape. 

To avoid such a systematic bias, studies of the DTD should 
strive to use consistent extinction corrections between low 
and high redshift. To aid this, we have provided our rates 
calculated under a variety of extinction assumptions. Con- 
sitency will go a long way towards reducing potential er- 
rors in the DTD, even if the extinction distribution remains 
poorly known. However, in the long run one would like a bet- 
ter understanding of SN extinction distributions at both low 
and high redshift, particularly to avoid uncertainties due to a 
changing extinction distribution with redshift. The prospects 
for directly detecting the highly extincted SNe are not great: 
even with a deep IR search, most SNe with Ay > 2 will 
be missed at high redshift. The alternative is more pre- 
cise updated modeling of SN la extinction in the vein of 
Riello & Patat (2005) or Mannucci et al. (2007) that takes into 
account factors such as the evolution of extinction with delay 
time and our latest knowledge of the SN la DTD. 

5.3. Summary 

In this paper we have computed volumetric SN la rates 
based on 189 HST orbits. This large HST dataset adds statis- 
tics to the existing HST rate measurements, previously based 
only on the GOODS fields. Our results provide additional 
strong evidence that the SN la rate is > 0.6 x 10~ 4 bjj yr -1 
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Mpc~ 3 at z ~ 1. The availability of raw data from our ef- 
ficiency simulations makes it simple to combine this dataset 
with current and future HST datasets, such as the 902-orbit 
Cosmic Assembly Near-IR Deep Extragalactic Legacy Sur- 
vey (CANDELS; Grogin et al. 2011), for even greater statis- 
tical gains. 

We find that an important systematic uncertainty in our re- 
sult is the amount of host-galaxy dust assumed in our simula- 
tions. This illustrates the need to use caution in comparing SN 
rate results from different surveys, especially as statistical un- 
certainty decreases and systematic uncertainties become dom- 
inant. Consistent comparisons and updated extinction models 
can drastically reduce the dust systematic in future studies. 
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